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ABSTRACT 

We present multiwavelength observations of the ultraluminous blazar-type ra- 
dio loud quasar PKS 0528+134 in quiescence during the period July to December 
2009. Four Targct-of-Opportunity (ToO) observations with the XMM-Newton 
Satellite in the 0.2 - 10 keV range were supplemented with optical observations 
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at the MDM Observatory, radio and optical data from the GLAST- AGILE Sup- 
port Program (GASP) of the Whole Earth Blazar Telescope (WEBT) and the 
Very Long Baseline Array (VLB A), additional X-ray data from the Rossi X-ray 
Timing Explorer (RXTE; 2-10 keV) and from Suzaku (0.5 - 10 keV) as well 
as 7-ray data from the Fermi Large Area Telescope (LAT) in the 100 MeV - 
200 GeV range. In addition, publically available data from the SMARTS blazar 
monitoring program and the University of Arizona / Steward Observatory Fermi 
Support program were included in our analysis. 

We found no evidence of significant fiux or spectral variability in 7-rays and 
most radio bands. However, significant fiux variability on a time scale of several 
hours was found in the optical regime, accompanied by a weak trend of spectral 
softening with increasing flux. We suggest that this might be the signature 
of a contribution of unbeamed emission, possibly from the accretion disk, at 
the blue end of the optical spectrum. The optical flux is weakly polarized with 
rapid variations of the degree and direction of polarization, while the polarization 
of the 43 GHz radio core remains steady, perpendicular to the jet direction. 
Optical spectropolarimetry of the object in the quiescent state suggests a trend of 
increasing degree of polarization with increasing wavelength, providing additional 
evidence for an unpolarized emission component, possibly thermal emission from 
the accretion disk, contributing towards the blue end of the optical spectrum. 
Over an extended period of several months, PKS 0528-1-134 shows moderate 
(amplitude < 50 %) flux variability in the X-rays and most radio frequencies 
on ~ 1 - 2 week time scales. We constructed four spectral energy distributions 
(SEDs) corresponding to the times of the XMM-Newton observations. We find 
that even in the quiescent state, the bolometric luminosity of PKS 0528-1-134 is 
dominated by its 7-ray emission. 

A leptonic single-zone jet model produced acceptable fits to the SEDs with 
contributions to the high-energy emission from both synchrotron self-Compton 
radiation and Comptonization of direct accretion disk emission. Fit parameters 
close to equipartition between the energy densities of the magnetic field and 
the relativistic electron population were obtained. The moderate variability on 
long time scales, compared to expected radiative cooling time scales, implies the 
existence of on-going particle acceleration, while the observed optical polarization 
variability seems to point towards a turbulent acceleration process. Turbulent 
particle acceleration at stationary features along the jet therefore appears to be 
a viable possibility for the quiescent state of PKS 0528-1-134. 

Subject headings: galaxies: active — Flat Spectrum Radio Loud Quasars: indi- 
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vidual (PKS 0528+134) — radiation mechanisms: non-thermal 



Introduction 



Blazars (BL Lac objects and gamma-ray loud flat spectrum radio quasars [FSRQ]) 
are the most extreme type of active galactic nuclei (AGN). They were historically defined 
through extreme flux variability throughout the electromagnetic spectrum, and sometimes 
strong and variable linear polarization at radio and optical wavelengths. In the 1990s, 
observations by EGRET on board the Compton Gamma-Ray Observatory revealed large 7- 
ray fluxes (often dominating the bolometric luminosity of the source) from many blazars. 
The radio through optical emission from blazars is commonly interpreted as synchrotron 
emission from ultrarelativistic electrons in a relativistic plasma jet that is closely aligned 

o 

with our line of sight {9obs < 20 ). This assertion is supported by the superluminal mot ion 



that most blazars exhibit (e.g., IJorstad et al.ll200ll : iLister et al.ll2009t iPiner et al.l 120061 ) as 



well as by the observed luminosity and variability timescales observed in these objects. In 
extreme cases, variability time scales down to a few minutes have been found in the very 
high energy (VHE) 7-ray regime (e.g.. 



Albert et al.l 120071 : lAharonian et al.ll2007n 



Two competing classes of models are currently being considered for the origin of the 
high-energy (X-ray through 7-ray) emission from blazars. In leptonic models, hadrons (pri- 
marily protons) in the jet (if present in substantial numbers at all), are assumed not to be 
accelerated to ultrarelativistic energies. They do not exceed the threshold for photo-pion 
production processes on the low-frequency radiation field in the jet, and proton synchrotron 
radiation is assumed to be negligible. Therefore, in leptonic models, the high-energy ra- 
diative output is dominated by Compton scattering of low-frequency photons off relativistic 
electrons. In hadronic models, it is assumed that ultrarelativistic protons exist in sufficient 
number. In such a scenario, the protons will dominate the radiative output via proton syn- 
chrotron radiation and synchrotron and Compton emission from secondary particles. Those 
are produced in photo-pion production and subsequent pion and muon decay an d electro- 
magne tic cascade processes. For a recent review of blazar emission models, see, e.g.. lBottcher 



torn . 



The mechanism(s) of acceleration of particles to ultrarelativistic energies in blazar jets 
are currently very poorly underst ood. Particle acceleratio n may be related to relativistic 
shocks in an unsteady flow (e.g., iMarscher fc Gearl Il985l ). internal sho cks resulting frona 



the collision of relativistic plasma blobs ej e cted at different speeds 



(e.g., 



Mimica et al. 2004:|joshi fc Bottcherll2010l : iBottcher fc Dermerll2010 ). re-coUimation shocks 



e.g. 



Spada et al.l 1200 1 



Bromberg fc LevinsonI l2009l ) . or relativistic shear layers in radially stratified jets (e.g. 
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Stawarz fc Ostrowskill2002l : iRieger fc Duffyll2004j . l2006l ). to name just a few plausible scenar- 
ios. Signatures that reveal the nature of particle acceleration in blazar jets may be found 
both in spectral and variability features. The nature of the acceleration mechanism is re- 
flected in the shape of the produced particle spectra. Those, in turn, can be inferred from 
the shape of the non-therma l photon spectra, in particular in the synchrotron part of the 
SED (e.g., iFinke et al.l 120081 ). The dynamics and light travel time effects, in particular in 
shoc k acceleration scenarios, w ill leave distinct imprints in the observed variability features 



[e.g. 



Bottcher k Dermei]l2010[ ). 



Observational studies of blazars have so far mostly focused on bright, flaring states of 
blazars. This is the consequence of observational constraints which make detailed measure- 
ments of spectral and variability features in X-rays and 7-rays difficult in low flux states. 
However, blazars are known to spend most of the time in their quiescent state which has so 
far received very little attention and is therefore very poorly understood. EGRET detected 
7-ray blazars almost exclusively in flaring states, and the simultaneously operating X-ray 
telescopes (ROSAT, ASCA, RXTE) lacked the sensitivity to measure detailed X-ray spec- 
tral and variability properties of most blazars in their quiescent states. Therefore, even the 
question whether 7-ray emission persists at all in the quiescent states of blazars remained 
an open issue during the EGRET era. 

The observational situation has dramatically changed with the advent of the new gener- 
ation of X-ray observatories, in particular Chandra and XMM-Newton as well as the launch 
of the Fermi Gamma-Ray Space Telescope in June 2008. The Fermi Large Area Telescope 
(LAT, lAtwood et al.l 120091 ) is continuously monitoring the entire sky every 3 hours in the 
energy range 20 MeV - 300 GeV with about an order of magnitude superior sensitivity 
compared to that of EGRET. It routinely detects 7-ray emission from known blazars even 
in their quiescent states. A detailed study of the quiescent state of blazars may elucidate 
whether the quiescent jet flow is smooth, exhibiting little or no variability, or the quiescent 
emission consists of the superposition of a rapid succession of "mini-flares". In particular, 
a featureless light curve in all bands might indicate the persistence of particle acceleration 
mechanisms not related to impulsive (shock) events, and might point towards shear-flow 
acceleration in radially stratifled jets, or standing features such as re-collimation shocks in 
the quiescent states of blazars. 

This situation has motivated us to propose Target-of-Opportunity (ToO) observations 
with XMM-Newton in AO-8, triggered by an extended quiescent state of a known 7-ray 
bright blazar. We deflned a quiescent state of a 7-ray blazar by the object maintaining a 
> 100 MeV 7-ray flux lower than the lowest flux or upper limit ever determined by EGRET, 
over at least 2 weeks. The prominent high-redshift 7-ray bright FSRQ PKS 0528+134 
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fulfilled our pre-specified trigger criterion through its continued 7-ray quiescence for several 
months prior to September 2009. We therefore triggered our XMM-Newton ToO observations 
on PKS 0528+134. The observations consisted of four observations on 2009 September 8, 
10, 12, and 14. These were coordinated with ground-based radio and optical observations. 
Simultaneous 7-ray observations were provided by Fermi LAT. 

In the following section, we give an overview of the known properties of our target, 
PKS 0528-1-134. In ^ observations and data reduction procedures are described. In ^ we 
present the results of a flux and spectral variability analysis. The structure of the parsec 
scale jet of this source is presented in section ^ Results of our modeling of four simultaneous 
SEDs obtained during our campaign are discussed in §6] A discussion on the optical spectral 
variability and other relevant issues is presented in section ^ Finally, we summarize our 
results and draw conclusions in ^|H1 Throughout this paper, we refer to a spectral index a 
as the energy index such that oc z/~", corresponding to a photon index Fph = a + 1. We 
use a ACDM cosmology with Qm = 0.3, Qa = 0.7, and Hq = 70 km s~^ Mpc~^. In this 
cosmology, the luminosity distance of PKS 0528+134 is cIl = 16.2 Gpc. 



2. The Quasar PKS 0528+134 



The compact FSRQ PKS 0528+134 is or ie of the most lumin ous and most distant 7-ray 
blazars known, with a redshift of 2; = 2.07 (IHunter et al.lll993l ). In the high-energy 7-ray 
band (abov e 100 MeV), this so urce was first detected by EGRET during the period 1991 
April- June (IMattox et al.lll997l ). Besides EGRET, this source was also detected by the other 
two instruments onboard of CGR O: the Oriented Scintillation Spectrometer Experiment 
(OSSE) in the 0.05 - 1.0 MeV band (IMcNaron-Brown et al.lll995l ). and the Imaging Compton 
Telescope (COMPTEL) 0.75 - 10 MeV). During EGRET obs ervations (from 1991 - 2000) 
PKS 0528+134 showed intense variability (iDingus et al.1 Il996h exhibiting strong flares in 
1991, 1993, 1995, and 1996. Its highest 7-ray flux was detected in 1993 March, when it 
reached 10^^ JyHz (IMukherjee et al.l 119961 ). strongly dominating the bolometric luminosity 
of the sour ce. PKS 0528+134 is faint in the optical with a mean visual magnitude of 
irty = 19.5 (jWall fc Peacockl 119851 ). This is a consequence of the high Galactic extinction 



in the direction o f PKS 0528+134, estimated to be 2.782 < Ay < 5 flSchleeel et al.l Il998 



Zhang et al.lll994l ). The reason for this high extinction is that PKS 0528+134 with Galactic 
coordinates / = 191.37°, b = -11.01° (a = 5^30 "^56^41, S = +13°31^55.^15 J2000), is locate d 



behind the translucent molecular cloud B30 (ILiszt fc Wilsonlll993l : iHogerheijde et al.lll995l ). 
Hence, there have been relatively few optical observations of this source compared to other 
bands. 
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In the radio regime, PKS 0528+134 is regularly monitored by several programs at dif- 
ferent frequencies. The source show s pronounced rad i o fiux-density variability on timescales 



of several months to a few years jAUer et al.l llQSsI: iReich et 



Stevens et al. I ll994i IValtaoja fc Terasrantalll995l : lPohl et al 



19961: 



1993 



Izhang et al 



Peng et al.ll2001 



1994 



Bach et al. 



20071 ). A delay from high frequencies to low frequencies in radio bursts has been identified, 
and delays of a few months betwe en 7-ray flares and the corresponding radio bursts have 
been found (IMukherjee et al.lll996l ). Additio nally, from VLBI ob servations performed in the 
8.4 GHz band over a period of almost 8 years iBritzen et al.l ( 1l999l ) found that PKS 0528+134 
has a bent jet of length ^5-6 mas extending toward the northeast, in which rapid struc- 
tural changes take place. Increasing activity in the radio and 7-ray bands is associated 
with morphological changes in the radio structure o f the jet, and superl uminal motion with 
/^±,app ^ 30 is found in some of the jet components (jjorstad et al.ll2005l ). 



Prior to the time period of EGRET, observations at X-ray energies were carried out 
by the Einstein Observatory in 1980, but no high confiden ce values for the X-ra y flux and 
spectral index were found due to the low count statistics (IBregman et al.l 119851 ). Starting 
in March 1991, and later in September 1992, X-ray observations in the energy range 0.07 - 
2.48 keV, performed with the ROSAT Position Sensitive Proportional Counter (PSPC) were 



used to investigate the georn etry and physical environment of PKS 0528+134 (jZhang et al. 



I994J : IMukherjee et al.lll996l ). The continuum emission of this source in the medium - hard 
X-ray band (0.4-10 keV) was first measured using observati ons with the Advanced Satellite 
for Cosmology and Astrophysics (ASCA) in 1994 and 1995 ( ISambruna et al.lll997l ). Further 
X-ray observations of PKS 0528+134 have been carried out with the Rossi X-ray Timing 
Explorer (RXTE) during August and September 1996 and May 1999, as well as with Bep- 
poS AX in the 0.1-10 keV and 15-200 keV bands during 1997 February and M arch as part of a 



multiwavelength campaign involving EGRET and ground based telescopes (iGhisellini et al. 
1999h . 



SEDs of PKS 0528+134 collected during six year s of EGRET observation s were compiled 
and modeled with a one-zone leptonic jet model by IMukherjee et al.l (Il999[ ). It was found 
that during all EGRET detections of the source, the bolometric luminosity was dominated 
by its 7-ray output. While, due to often incomplete multiwavelength coverage, the modeling 
results of that paper were subject to a large degree of freedom and uncertainty, the observed 
epoch-to-epoch variability of PKS 0528+134 was found to be consistent with a correlation 
between the 7-ray flux and the bulk Lorentz factor of the emission region along the jet. 

Given the extreme properties of PKS 0528+134, this blazar has been the target of many 
observations at different wavelengths. However, as for almost all blazars (see coordinated 
multiwavelength campaigns have targeted flaring states. The quiescent-state SEDs and 
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spectral variability patterns of blazars in general, and of PKS 0528+134 in particular, are 
poorly understood. Given its extended 7-ray quiescence throughout 2009, as revealed by 
Fermi LAT minotorin^, PKS 0528+134 was therefore found to be an appealing target for 
our pre-approved XMM-Newton Cycle 8 ToO observations and multiwavelength campaign. 



3. Observations and Data Reduction 

The blazar PKS 0528+134 was the target of intensive, simultaneous and quasi- simultaneous 
observations at optical (MDM, GASP), radio (GASP), X-ray {XMM-Newton, RXTE, Suzaku), 
and 7-ray [Fermi LAT) frequencies during the period September 8 to 18, 2009. In addition, 
a more extended period of time, throughout July - December 2009, was covered by less 
intensive radio, optical. X-ray, and 7-ray monitoring for longer-term (weeks - months) vari- 
ability studies. In addition to these previously unpublished data, we included publically 
available photometric monitoring data from the Small and Moderate Aperture Research 
Telescope System (SMARTS)! at the Cerro Tololo Interamerican Observatory, in our data 
collection. Specifically, we included B, V, and R-band data from the Yale Fermi/SMARTS 
project, covering the entire, extended campaign period. We also included publically available 
polarimetry and spectroscopy data from the University of Arizona - Steward Observatory 
Fermi support progranjfl in our analysis. 



3.1. Optical Observations 

3.1.1. Optical Photometry 

Most of our optical (BVR) observations were carried out using the 1.3-m McGraw-Hill 
telescope of the MDM Observatory on the south-west ridge of Kitt Peak, Arizona. The 
telescope is equipped with a 1024x1024 pixels CCD camera and standard Johnson- Cousins 
UBVRI filters. Every night during the period September 9 - 19, 2009, sequences of science 
frames on PKS 0528+134 were taken in the B-V-R filters with exposure times of 360, 180, 
and 180 seconds, respectively, with slight variations depending on atmospheric conditions. 
All frames were bias-subtracted and flat-field corrected using standard routines in IRAF. In 



^see http://fer111i.gsfc.nasa.g0v/ssc/d.ata/access/lat/mslJ.c/ 



http : //www . astro . yale . edu/ smarts/ 



http : // j ames . as . ar izona . edu/~psmith/Fermi/ 
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most cases, the signal-to-noise ratio (SNR) of PKS 0528-1-134 on individual, reduced frames 
was too low to allow for a high-precision magnitude determination. Due to the high redshift 
oi z = 2.07 and the expectation that the time scale (in the cosmological rest frame of the 
blazar) of brightness variations in FSRQs is > 1 hr, we typically added reduced frames taken 
within intervals no longer than three hours in order to increase the SNR. 

In order to construct light curves of this quasar we applied re lative photornetry t o the 
frames using the three comparison stars in the sky chart provided by lRaiteri et al.l (119981 ). We 
extracted instrumental magnitudes of the three comparison stars and PKS 0528+134 using 
the phot routine within the DAOPHOT package of IRAF. The light curves thus constructed 
are analyzed in section §H 

Additional BVRI photometric and R-band polarimetric observations of PKS 0528+134 
were performed at the 1.8 m Perkins telescope of Lowell Observatory (Flagstaff, AZ) during 
a 2-week campaign from October 15 to October 28 2009 with the PRISM camereo. Stan- 
dard bias subtraction and flat field correction were carried out for each frame. For these 
images we also performed differential photometry with an aperture o f radiu s 6" using BVR 
measurements for comparison stars 3, 2, and 1 from iRaiteri et al.l (119981 ). We obtained 
I-band measurements for the same comparison stars (14.475 ± 0.010, 13.239 ± 0.010, and 
12.220 ± 0.010 ma g, respectively) using comparison stars from the field of PKS 0735+138 
(jSmith et al.l Il985l ) observed just after PKS 0528+134 on 7 nights. The PRISM camera 
possesses a polarimeter with a rotating half wave plate that we employed for R-band po- 
la rimetry. Most details of the polarization observations and data reduction can be found 



m 



Jorstad et al.l (120101 ). The polarization data were corrected for the statistical bias associ- 



ated with the fact that th e degree of polarization is restricted to being a positive quantity 
JWardle fc Kronberel 119741 ). 



We also obtained optical data taken using the polarimetric imaging system Polima 
attached to the 84 cm telescope at San Pedro Martir, Baja California, Mexico as part of a 
long-term polarimetric monitoring program of a sample of 35 blazars during 6 photometric 
nights in October and November 2009. Due to the optical design of Polima, 4 images with 
different position angles of the polarimeter to measure the Stokes parameters were taken 
through an R-band filter. The data reduction for these frames (correction for bias and pixel- 
to-pixel variations across the CCD) as well as the photometry has been carried out on each 
of the individual images using a dedicated pipeline developed by D. Hiriart. The fluxes 
of PKS 0528+134 and several stars in the field of view for each PA have been combined. 
Flux calibration was finally done via the comparison stars 2 and 3 in the field of view from 



"^http: / /www. bu.edu/prism / 
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Raiteri et al.l fllQQSh 



3.1.2. Optical Spectroscopy 

In addition to the photometric observations described in the previous sub-section, PKS 0528+134 
is also regularly monitored with spectroscopic observations through the University of Arizona 
- Steward Observatory Fermi Support program. As a representative example, we show in 
Figured], the optical spectrum of October 19, 2009, which is within the extended period of 
our multiwavelength campaign. The source was clearly still in the quiescent state targeted 
in this work. 

The spectrum exhibits two distinct emission lines: GUI] A1909, red-shifted to A = 
5860 A, and CIV A1550, red-shifted to A = 4740 A. The total measured flux in the lines 
corresponds to Fcm] = 1.8 x 10~^^ erg cm~^ s~^ and Fqiv = 4.3 x 10"^^ erg cm~^ s~^. In 
order to evaluate the de-absorbed fluxe s, we evaluate the extinction coefficients Ax with the 
extinction law of ICardelli et al.l (Il989[ ) and Ay = 2.78, yielding ^4740 = 3.35 and A^sqq = 
2.59. This yields intrinsic luminosities in the two lines of Lcm] ~ 5.8 x 10^^ erg s~^ and 
Lciv ~ 2.5 x 10^^ erg s-\ 

The continuum redward of ~ 5000 A can be well fit with a power-law Fx oc A^'^, which 
corresponds to a steep spectrum in frequency space, Fi, oc z/~^'^. 



Optical Polarimetry 



R-band photo-polarimetric observations of PKS 0528+134 were also acquired with the 2.2 m 
telescope of the Calar Alto Observatory in Ameri'a, Sp ain, as part of the M APCAtB program. 
The d ata were reduced and calibrated as described in Agudo et al. (j2011 ) and Jorstad et al. 

feoioh . 



St. Petersburg observations were performed at the 70-cm telescope of Crimean Obser- 
vatory using ST-7 XME photometer-polarimeter. The standard procedure, including bias 
and dark s ubtraction, flat-field correction and calibration relative to comparison stars 1, 2 
and 3 from iRaiteri et al.l (119981 ) was applied. 



Additional polarimetry and spectropolarimetry data were included from the University 
of Arizona - Steward Observatory Fermi Support program. Those data included two spec- 



^ Monitoring AGN with Polarimetry at the Calar Alto Telescopes: 
http: // www, iaa. es / ^iagudo /research/M AP CAT] 
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C IV ^ 1550 




X [Angstrom] 



Fig. 1.— Optical spectrum of PKS 0528+134 on December 19, 2009, taken at the Steward 
Observatory. The two prominent emission hnes of GUI] and CIV are labeled and correspond 
to luminosities of Lcm] ~ 4.2 x 10^^ erg s~^ and Lciv ~ 3.2 x 10^^ erg s~^. The continuum 
red- ward of ~ 5000 A can be fit with a power-law Fx oc A^-^. 
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Fig. 2. — Wavelength dependence of the degree of polarization of the optical emission of 
PKS 0528+134, from two observations in 2008 December and 2009 January. The plot is 
highly suggestive of a systematic trend of increasing polarization with increasing wavelength. 
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tropolarimetric observations of PKS 0528+134 in 2008 December and 2009 January, which 
showed a rather high (P > 10 %) degree of polarization, allowing for a meaningful spectrap- 
olarimetric analysis. At those times, the object was in a similar quiescent state as during 
our campaign. The contributions of the CIII] and CIV lines (see previous sub-section) were 
assumed to be unpolarized, and their fluxes subtracted from the respective energy bins. 
Figure [2] shows the degree of polarization as a function of wavelength. The plot is highly 
suggestive of a systematic increase of the polarization towards longer frequencies. This may 
be interpreted as an increasing contribution of synchrotron emission towards the red end of 
the optical spectrum. 



3.2. X-Ray Observations 



3.2.1. XMM-Newton 



PKS 0528+134 was observed by the XMM-Newton Observatory (jjansen et al.l l200l[ ) 
between September 8th and 14th 2009 on four consecutive revolutions. Table [1] summarizes 
these observations. Here, we focus on the data taken with the EPIC detector, covering the 
energy band between 0.2 and 10 keV. EPIC observations were taken in full frame mode and 
thin optical filter (except the pn exposure in 0600121501 that was tak en on small window 
mode). XMM-Newton also has an Optical/UV Monitor Telescope (OM) jMason et al.lboOlj ). 
a small 30 cm telescope co-aligned with the main XMM-Newton X-ray telescopes. However, 
due to the optical faintness of PKS 0528+134 (> 19*^^ mag), OM observations did not return 
useful data. 



Obs. ID 


Date 
yy/mm/dd 


Exp. (pn) 
(ksGc) 


CR (pn) 
(cts/sec) 


Exp. (MOSl) 
(ksec) 


CR (MOSl) 
(cts/sec) 


Exp. (MOS2) 
(ksec) 


CR (MOS2) 
(cts / sec) 


0600121401 
0600121501 
0600121601 
0600121701 


2009-09-08@05:55UT 
2009-09-10a06:05UT 
2009-09-12903:01UT 
2009-09-14ia02:52UT 


9.94(14.28) 
19.89(20.13) 
21.11(21.96) 


0.2675±0.0052 
0.2194±0.0033 
0.2616±0.0035 


15.2(22.1) 
23.63(28.46) 

26.5(26.5) 
26.68(35.61) 


0.0897±0.0024 
0.0848±0.0019 
0.0808±0.0017 
0.0754±0.0017 


15.17(21.49) 
24.38(28.48) 

26.6(26.6) 
29.72(36.12) 


0.0943±0.0025 
0.0897±0.0019 
0.0908±0.0018 
0.0877±0.0017 



Table 1: Summary of PKS 0528+134 XMM-Newton EPIC observations. Exp. = Livetime 
after (before) correction due to periods of higher background activity. CR = Count rates, 
given for the 0.2-10 keV energy range. The quadrant of the pn chip containing the source of 
the pn exposure in observation 0600121701 was lost and hence there is no source information. 



Data Reduction 



The data have been analysed using SASv9.0 (iGabriel et al.ll2004j ) and corresponding calibra- 
tion files. Event and source lists were obtained for the EPIC detector, following the standard 
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SAS data reduction procedures. Several filtering criteria have been applied. The event list 
has been filtered for time periods of high background activity following the standard proce- 
dure of removing those time periods with background event rates at energies E > 10 keV 
higher than 1.0 cts/sec and 0.8 cts/sec for EPIC-pn and EPIC-MOS, respectively. The time 
losses due to the removal of these periods are up to 30 % depending on the observation 
and the instrument (three out of the four observations were performed at the end of their 
respective revolutions with the consequent increase in radiation levels towards the end of the 
observations). Table [T] shows the livetimes for each observation and instrument before and 
after this correction has been performed. The data was also filtered to include only single and 
double (PATTERN<4) pattern events for EPIC-pn and single to quadruple (PATTERN<12) 
for EPIC-MOS as well as those with quality FLAG=0. The filtered event lists were used 
to generate light curves and spectral products. The source region was considered to be a 
circular region centered around the source. Annular background regions were chosen to be 
centered around the source with radii 60" < R < 80". The source extraction regions were 
then optimized based on S/N for the given background region, which yielded typical values 
of ~ 40" in radius for the source region. 



EPIC Spectral Analysis 

Time-averaged spectra have been obtained for each individual observation. The spectra were 
re-binned in order not to oversample the intrinsic energy resolution of the EPIC cameras 
by a factor larger than 3, while making sure that each spectral channel contains at least 25 
background-subtracted counts. Both conditions allow the use of the quality-of-fit esti- 
mator to find the best fit model. Fits were performed in the 0.2 - 10 keV energy range. 
The spectra from all three EPIC cameras have been used simultaneously during the fitting 
procedure. The systematic difference between the EPIC cameras is belo w ~ 5 % in nor mal- 



ization. For the spectral analysis and fitting procedure, XSPEC vl2.4 (jArnaudlll996[ ) was 
used. 

The spectral model used to fit the data is composed of a power law convolved with a 
combination of Galactic column density of absorbing material traced by 21-cm emission from 
atomic hydrogen (NH,2icm), plus an additional absorbing column density (NH,moi), primarily 
due to heavy elements in the intervening molecular cloud B30, which are not properly traced 
by the hydrogen column density. The Galactic hydrogen column density is kept fixed during 
the fitting procedure. The spectral fitting model takes the following form for the differential 
photon flux ^{E): 
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where cr{E) is the photoelectric absorption cross-section, with abundances after I Anders fc Grevesse 



( 119891 ). and Fph the power-law index with normalization A^. Errors in the relevant parameters 



are given at the 90 % confidence level (CL) for any given parameter. 

Table [2] shows the best fit values for the model considered for all four observations. 
Figure [3] shows an example of the spectra as determined for one of the observations (composed 
of EPIC-pn and the two EPIC-MOS spectra) with the best fit model and deviation. 
Similar spectra are derived for the other observations. 



3.2.2. RXTE 

In addition to our four ToO XMM- Newton pointings, we monitored the 2.4 - 10 keV 
X-ray fiux of PKS 0528+134 with the Rossi X-ray Timing Explorer (RXTE) Proportional 
Counter Array (PCA), with exposure times of 900 - 3400 s for individual observations. The 
X-ray fiux measurement entailed the subtraction of an X-ray background model (faint source 
model, version 20051128) from the raw spectrum using the standard X-ray data analysis 
software packages FTOOLS and XANADU. We used the program XSPEC to fit the residual 
photon spectrum with a power-law model plus photoelectric absorption along the line of 
sight. 



3.2.3. Suzaku 

PKS 0528-1-134 was also observed with the Suzaku satellite as a part of multi-band 
observations conducted in September 2008, about one year prior to the campaign covered 
in this paper. However, both the X-ray fiux and 7-ray fiux levels were comparable to those 
measured during the September 2009 campaign, indicating a quiescent state. We therefore 
include the Suzaku data here for comparison of the spectral and short-term variability prop- 
erties with those measured by XMM-Newton one year later, in a similar quiescent state. The 
Suzaku satellite features instruments sensitive in the soft X-ray band. 

Suzaku observations of PKS 0528+134 started on 2008 September 27, 02:38 UT, and 
lasted until 2008 October 2, 16:12 UT. Since the source was not detected in the Hard 
X-Ray Detector (HXD) data, we considered only the X-Ray Imaging Spectrometer (XIS) 
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PKS0528+134 EPIC Spectrum (0.2-10 keV) 



0.01 




10-3 



1 2 
Energy (keV) 

Fig. 3.— Example of EPIC (pn, MOSl and M0S2) spectrum for the observation id 
0600121601 where EPIC data is available for all three cameras. The fit has been done 
in the 0.2-10 keV energy range with channel grouping > 3 channels and > 25 background 
subtracted counts. 



-17- 



instruments. The observation conditions were nominal, although the XISl data suffered from 
unusually high and variable background, resulting in a total apparent background count rate 
ranging from 1 to 3 counts s~^ over the entire chip. Nonetheless, since the background- 
subtracted spectrum determined from the XISl data was entirely consistent with that from 
XISO and XIS3, we included the properly background-subtracted XISl data in the spectral 
fitting (see below). 

The total exposure time yielding good data accumulated in the pointing was 203 ks. 
We used the standard ftools data reduction package, provided by the Suzaku Science 
Operations Center, with the calibration files included in the CALDB ver. 4.3.1. The net 
count rates were, respectively, 0.10, 0.13, and 0.12 count s~^ for XISO, XISl, and XIS3. For 
the analysis of spectra and light curves, we extracted the counts from a region corresponding 
to a circle with 260 arc sec radius. We used a region of a comparable size from the same chip 
to extract the background counts. We plot the resulting total light curve (not background- 
subtracted) from XIS3, binned in 5400 sec bins, in Figure HJ and note that the background 
was steady at the level of ~ 0.015 count s~^. The data indicate no significant rapid (a day or 
less) variability during the Suzaku observation, but show a long-term trend, with a ~ 20 % 
decrease in flux over the ~ 5.3 day long duration of the observation. A similar secular flux 
decrease is seen in the two other XIS detectors. 



For the spectral analysis, we we used the XSPEC spectral analysis software f lArnaud 



19961 ). For the spectral fitting, we used the standard redistribution files and mirror effec- 
tive areas generated with S'nzaA:u-specific tools. We included the counts corresponding to 
the energy range of 0.5 - 10.0 keV in our spectral fits. We used all three XIS detectors 
simultaneously, but allowed for a small (a few %) variation of the relative normalizations. 
The source spectrum was model ed as an absorbed power law, w ith the cross-sections and 



elemental abundances as given in lMorrison fc McCammoru (119831 ). Other absorption models 



give similar results. The best-fit absorbing column was 4.6 ±0.2 x 10^^ cm~^, and the photon 
index was 1.54±0.03, with the acceptable best-fit of 3407 for 3512 d.o.f. The Suzaku spec- 
trum is shown in Figured The observed model 2-10 keV flux is 2.8 x 10^^^ erg cm~^ s~^, 
with the statistical error of 5 %, which is probably smaller than the systematic error resulting 
from the calibration uncertainty of the Suzaku instruments. 

The fitted values of the total absorbing column are only marginally consistent between 
the Suzaku and XMM-Newton data sets. This might be in part related to the differences 
in calibration of the two instruments. An additional source for the discrepancy might be 
the fact that (1) we are using Galactic elemental abundances, which might or might not 
be appropriate for the line-of-sight molecular cloud, and (2) the bandpass of XMM-Newton 
extends down to 0.3 keV, while for Suzaku it is 0.5 keV. If the elemental abundances in 
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PKS 0528+134 Suzaku XIS3 observation 
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Fig. 4.— Suzaku X-ray light curve of PKS 0528+134 from the XIS3 data of 2008 September 
27 - October 2. t ^ corresponds to MJD 54736.1771. 
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Suzaku XISO, XISl and XIS3 spectrum of PKS 0528+134, Sep 2008 




1 2 5 

Energy (keV) 



Fig. 5. — 0.5 - 10 keV X-ray spectrum of PKS 0528+134 from the Suzaku observations of 
2008 September 27 - October 2. 
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the intervening absorber of low-z elements are not the same in the molecular cloud as the 
Galactic abundances assumed in the spectral fit, the precise value of the best-fit might 
be dependent on the bandpass of the instrument and even the details of the effective area at 
the energies where the absorption edges due to those elements play an important role. 

While the X-ray spectral index measured by Suzaku in 2008 September/October is 
consistent with the value measured one year later by XMM-Newton (see §3.2. H table [2]), the 
Suzaku 2-10 keV flux is about a factor 2.4 higher than the XMM-Newton value. This is in 
agreement with the ~ factor 2-3 variability measured by RXTE on ~ weekly time scales 
(see gH). 



3.3. 7-Ray Observations 

Gamma-ray observations of PKS 0528+134 were performed by the Fermz-LAT. This is 
a pair-conversion 7-ray telescope sensitive to photon energies greater than 20 MeV. In its 
nomin al scanning mode, it surveys the whole sky every 3 hours with a field of view of about 



2.4 sr (lAtwood et al.ll2009l ). The data presented in this paper (restricted to the 100 MeV 
- 200 GeV range) were collected from MJD 54983 (2009 June 1st) to MJD 55193 (2009 
December 28). For this analysis, the Diffuse event class was used. This is the optimized 
class for point source analysis with minimal residual contamination from charged-particle 
backgrounds. To minimize systematics, only photons with energies greater than 100 MeV 
were considered in this analysis. In order to limit the contamination from atmospheric 7-rays 
produced by interactions of cosmic rays with the upper atmosphere of the Earth, only events 
with zenith angle < 105° were selected. In addition, time intervals during which the rocking 
angle was larger than 52° have been excluded from the analysis, because the bright limb of 
the Earth enters the field of view. The analysis was performed with the standard analysis 
tool gtlike, part of the Fermi-LAT Science Tools software package (version v9rl5p5)|^. The 
P6_V3_DIFFUSE set of instrument response functions was applied. 

Photons were selected in a circular Region Of Interest (ROI) 7° in radius, centered at 
the position of PKS 0528+134. The isotropic background, including the sum of residual 
instrumental background and extragalactic diffuse 7-ray background, was modeled by fitting 
this component at high galactic latitude (file provided with Science Tools). The Galactic 
and Isotropic diffuse emission models version "gll_iem_v02.fit" and "isotropicJem_v02.txt" 



^For a documentation of the Science Tools see 



http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/ 
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were usedEI All point sources in the first Fermi-LAT catalog ( Abdo et al. 2010c ) lying 
within the ROI and a surrounding 10°-wide annulus were considered in the fit and modeled 
with single power-law distributions of the form F{E) = Nq{E / Eq)~^ . Their flux was kept 
free whereas their spectral index value was frozen to the value listed in the IFGL catalog, 
except for PKS 0528+134 whose index was kept free. Due to limited statistics, the gamma- 
ray spectrum of PKS 0528+134 was modeled with a power-law distribution in the present 
analysis, although the spectrum measured over 11 months covered by the IFGL catalog 
exhibits distinct curvature, characterized by a curvature index of 8.14, corresp onding to a 



4 % p robabihty that the spectrum is adequately represented by a power law f lAbdo et al. 



2010aJ). The highest energy photon attributed to PKS 0528+134 has an energy of 8.6 GeV. 
The estimated systematic uncertainty on the flux is 10 % at 100 MeV, 5 % at 500 MeV and 
20 % at 10 GeV. 



3.4. Radio Observations 

We obtained radio observations from the GLAST- AGILE Support Program (GASP) of 
the Whole Earth Blazar Telescope (WEBT). In table |3l we list the frequencies in which we 
collected data and the observatories participating in the campaign. 

Centimeter-band observations were obtained with the University of Michigan 26-meter 
prime focus paraboloid equipped with radiometers operating at central frequencies of 4.8, 
8.0, and 14.5 GHz. Observations at all three frequencies utilized rotating polarimeter sys- 
tems permitting both total flux density and linear polarization to be measured. A typical 
measurement consisted of 8 to 16 individual pointings over a 20 - 40 minute time period. 
Frequent drift scans were made across stronger sources to verify the telescope pointing cor- 
rection curves. Observations of program sources were intermixed with observations of a grid 



flux scale was based on observations of Cassiopeia A fse e. e.g.. iBaars et a. 



the calibration and analysis techniques are described in lAUer et al.l ( 119851 ) 



^1 ^1g-1 


Baars et a 


1977) 


AUer et al. 


(1985) 





The 37 GHz observations were made with the 13.7 m diameter Metsahovi radio telescope. 
The detection limit of the telescope at 37 GHz is of the order of 0.2 Jy under optimal 
conditions. Data points with a SNR < 4 are handled as non-detections. The flux density 
scale is set by observations of DR 21. Sources NGC 7027, 3C 274 and 3C 84 are used as 
secondary calibra t ors. A detailed description of the data reduction and analysis is given in 
Terasranta et al.l ( 119981 ). The error estimate in the flux density includes the contribution 



''For details on the background model see |http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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EPIC Power Law 
PHABS+PHABS 



Obs ID 



Nn.mol 



r 



ph -T 0.2-2keV 

10^^^ ergs cm~^ s^^ 



F2-10keV X /d.o.f. 

10~^^ ergs cm~^ s~^ 



0600121401 
0600121501 
0600121601 
0600121701 



1 5+°-°^ 
1 7+0-02 

'-'•J- ' -0.02 
1 2+0-64 



1.57: 

i.6o: 
i.6i: 
1.58: 



-0.07 
-0.07 
-0.05 
-0.07 
-0.05 
-0.02 
-0.07 
-0.10 



0.28: 
0.24: 
0.24: 
0.26: 



-0.03 
-0.04 
-0.02 
-0.03 
-0.02 
-0.02 
-0.03 
-0.02 



(0.88) 
(0.76) 
(0.81) 
(0.69) 



1.35 
1.10 



-1-0.11 
-0.11 
4-0.07 



-0.09 
-'-■-'■^-0.06 
^•^^-0.11 



(1.39) 0.93 (146.2/157) 

(1.14) 1.02 (243.3/238) 

(1.18) 1.05 (275.2/261) 

(1.22) 0.99 (135.6/137) 



Table 2: Fit performed in the 0.2 - 10 keV energy range. All 3 EPIC instruments have 
been used simultaneously in the fit. The Galactic absorption traced by 21-cm emission from 
neutral hydrogen, NH,2icm has been fixed to a value of 0.24 ■ 10^^ cm~^ according to the 



LAB Survey of Galactic Hl lKalberla et al.l (120051 ). Fluxes given are not absorption corrected 
(in parenthesis the de-absorbed flux is given). Errors indicate the 90 % CL. NOTE: For 
observation id 0600121701 only MOSl and M0S2 are available. In this case, the fit has been 
performed in the 0.3 - 10 keV energy range and the fluxes given are between 0.3 - 2.0 and 
2.0 - 10 keV. 



Table 3. Frequencies and observatories used to collect radio data 



Frequency Observatory 

05 GHz UMRAO, MEDICINA 

08 GHz UMRAO, MEDICINA 

14 GHz UMRAO 

22 GHz MEDICINA 

37 GHz Metshovi Radio Observatory (KURP-GIX) 

43 GHz Noto 

230 GHz MAUNA KEA (SMA) 

345 GHz MAUNA KEA (SMA) 
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from the measurement rms and the uncertainty of the absolute cahbration. 



The 43 GHz Noto observations have been performed using the On The Fly (OTF) scan 
technique (scan duration about 20 s) and the telescope gain (K/Jy) was determined as a 
function of the elevation using NGC 7027 as primary calibrator. To improve the signal to 
noise ratio many scans have been acquired and then averaged for a total integration time 
of about 15 minutes. The antenna temperature has been estimated by a gaussian fi t of th e 
average scan. A more detailed description for the 43 GHz data is given by lLeto et al.l ( l2009l ). 
Radio observations at the Medic ina radio ob s ervat ory were performed at 5, 8, and 22 GHz, 
and were analyzed as detailed in lBach et al.l (120071 ). 



Data at 230 GHz and 345 GHz were collected at the Submillimeter Array (SMA) on 
Mauna Kea and reduced using the MIR data reduction software. SMA flux density measure- 
ments are produced from a mixture of dedicated flux calibration/monitoring observations 
and data from science projects that may utilize a quasar as a calibration standard (typically 
for gain calibration of the interferometer). Raw visibility data are calibrated to correct for 
atmospheric absorption and instrumental gain variations, and then referenced to observa- 
tions of standard flux calibration sources, generally planets and/or moon s . Mor e details on 
the SMA flux density monitoring program can be found in iGurwell et al.l (120071 ). 



4. Variability Analysis 

One of the main goals of this multiwavelength campaign was the search for flux and 
spectral variability of PKS 0528+134 in its quiescent state. We first describe our results on 
flux variability in §4.11 and then turn to the investigation of spectral variability in §4.2[ In 
§4.31 we discuss variability of the optical polarization. 



4.1. Flux Variability 

The optical light curves of PKS 0528+134 obtained with the 1.3-m McGraw-Hill tele- 
scope of the MDM Observatory during the core of our multiwavelength campaign (September 
9 - September 19, 2009), along with the publically available SMARTS data from the same 
period, are plotted in figure [61 An increase in the brightness of the source of ~ 0.7 mag- 
nitudes in the R and V filters around JD 2455086 (September 11) is evident. In order to 
quantify the variability, we computed the reduced for a fit to a constant flux. Variabil- 
ity is evident in all three bands with = 15.9, 5.2, and 26.2, respectively, for the R, V, 
and B bands. Due to the limited observability period of PKS 0528+134 in any given night 
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(typically < 4 hr), the light curve during this mini-flare is clearly undersampled. Due to the 
sparse sampling of the optical observations, we can only constrain the variability time scale 
to tX* < 1 day. 

In the same way, we analyzed the X-ray variability on time scales of ~ 1-2 day 
between our four XMM-Newton observations as well on time scales of a few hours, within 
the individual observations. The XMM-Newton light curve is plotted in the bottom panel of 
figure O No flux variability is evident. For the 0.2 - 10 keV X-ray flux, a fit to a constant 
results in a of 0.91, thus confirming the absence of significant variability. For the same 
XMM-Newton observations, we also analyze the intraday variability. For most observations 
the for a fit to a constant flux is less than 1. In figure [3, we show the light curve from 
the XMM-Newton MOSl data of September 10, 2009, as an example. For this particular 
observation, ^ 0.81 is found. 

A multiwavelength variability analysis for a more extended period of time was made 
possible by including RXTE monitoring data. In figure [H we show the light curves of 
PKS 0528-1-134 in radio (5 and 37 GHz), optical. X-rays (RXTE), and 7-rays (bottom to 
top). The vertical shaded band highlights the interval of the core campaign where the 
most intense optical and X-ray observations (including XMM-Newton) were performed. The 
RXTE X-ray data from the period July 14 to December 2, 2009, represent the most extended 
coverage. In this interval, PKS 0528-1-134 shows significant variability corresponding to a 
reduced of 3.83 for a fit to a constant flux. The RXTE light curve indicates variability 
with flux changes of |AF/F| ~ 50 % on a characteristic time scale of t:^^, ~ 1 week. 

We note here that the fluxes resulting from our RXTE analysis are systematically higher 
than those measured by XMM-Newton during the same period. We have carefully double- 
and cross-checked both analyses and confirmed this discrepancy. The ROSAT All Sky Survey 
Catalogue does not list any known source in the field of view of RXTE which may be 
responsible for the flux discrepancy. Furthermore, three RXTE observations were carried 
out within a few hours of one of the XMM-Newton observations. Given this short time 
period and the systematic offset between the two instruments, rapid variability appears 
unlikely to cause the flux discrepancy. The factor of 4 discrepancy seems also too large to 
be solely due to calibration uncertainties. A plausible explanation for the systematic offset 
could lie in an uncertain background model for the RXTE analysis in the region around 
PKS 0528-1-134, near the Galactic anti-center. While this would cause a constant offset 
of the RXTE flux levels, it will not affect the variability. We are therefore confident that 
the RXTE flux variability analysis presented here is robust. For the analysis of the spectral 
energy distributions in ^ the flux and spectral information from XMM-Newton will be used. 

The optical light curve (second panel from the bottom of Figure [8]), including data 
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Fig. 6.— MDM and SMARTS optical light curves (RVB) of PKS 0528+134 between Septem- 
ber 9 and 19, 2009. The bottom panel shows the X-ray light curve corresponding to the four 
XMM-Newton observations. 



-26- 



0.20 




10000 15000 
Time [s] 



25000 



Fig. 7.— XMM-Newton MOSl light curve of our ToO observation on September 10, 2009. 
Each point corresponds to a bin of 500 seconds. No evidence for intra-day variabihty is 
found. 
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from the MDM, SMARTS, Perkins, San Pedro Martir, Calar Alto, and Crimean observato- 
ries, covers the interval from September 9 to November 19, 2009. This is the band where 
PKS 0528+134 shows the most significant variations in its flux density, with flaring episodes 
exhibiting brightness changes of AR < 1™^^ on a time scale of < 1 day. 

As shown in figure [8] (top panel), the Fermi 7-ray flux does not present significant 
variability during the extended campaign period (2009 July - December). The best fit found 
for the flux (E > 10 MeV) is (0.12 ±0 02) x 10"^ ph cm'^ s"^ with a test statistic TS = 68. 



The Test Statistic (iMattox et al.lll996l ) is defined as twice the difference in log(likelihood) 
obtained by including the source of interest and omitting it in the source model used in 
the gtlike analysis. This flux i s slightly lower tha n the mean flux observed over the first 



eleven months of Fermi data ( lAbdo et al.l l2010al ). which corresponds to (0.17 ± 0.01) x 



10 ^ ph cm ^ s ^ as represented by the dashed horizontal line in the same panel. 

Among the frequencies monitored in the radio regime, the 5 GHz and 37 GHz light 
curves had the most extended coverage in time. In the bottom panel of figure El we show 
both light curves. Only moderate variability with amplitudes of |AF/F| < 20 % is observed 
at these frequencies. In particular, a decreasing flux tendency is found at the end of the 
37 GHz light curve. A fit to a constant flux results in xl = 3.92 and 2.61 for the 5 GHz and 
37 GHz light curves, respectively. 

A more detailed analysis of the radio regime, including the light curves in all the mon- 
itored radio frequencies, is presented in flgure M The radio light curves show moderate 
variability in general. The 8 GHz and 14 GHz bands present signiflcant flux variations with 
xl = 32.96 and 10.17, respectively, for a flt to a constant flux. The results of a similar 
analysis for all radio frequencies are summarized in Table HI Variability appears to occur on 
time scales of t^lr'° ~ 1 ^ 2 weeks, comparable to the X-ray variability time scale measured 
by RXTE (see above). 

In summary, the comparative flux variability analysis of PKS 0528+134 at different 



Table 4. Light curve variability analysis in radio frequencies: xl fo^^ fits to a constant flux. 



Frequency xi 



5 GHz 3.92 

8 GHz 32.96 

14 GHz 10.17 

37 GHz 2.61 

230 GHz 0.62 



- 28 - 



0.5 
0.4 
i 0.3 
0.2 



2 0.1 1- 



6e-12 
g) 4e-12 
^ 2e-12|- ♦ 



= — ^ — \ — ' — \ — ^yn — ' — I ■ I ~ 

L , i i xi • 2.4- lOkeV 



Ob 



I 18.5 



D 1Q - 

S 19.5 

3.5 
^3.0 
^ 2.51- 
g 2.0 

Ei: 1-5 

1.0 



Fermi LAT (E > 100 MeV) 
1 1 months mean 



1 



\ =1 

♦ -_ 



X 



San Pedro 

MDM 

Perkins 

Calar Alto 

Crimean 

SMARTS 



X 



I I ! I I I I I I I I I 



— \ — ^ — r — ^ — \ — ^ — \ — ' — \ — ^ — \ — ^ — r 



I. 



1 



J ^ \ ^ L_ 



• 5 GHz 

♦ 37 GHz 



1 



1 



5020 5040 5060 



5080 5100 5120 
JD - 2,450,000 



5140 5160 



Fig. 8. — Light curves of PKS 0528+134. From top to bottom: a) the Fermi 7-ray flux 
in 15-day integration bins, b) X-ray (RXTE) light curve, c) Optical (R-magnitude) light 
curve including data from the MDM, SMARTS, Perkins, San Pedro Martir, Calar Alto, and 
Crimean observatories, d) Radio light curves at 5 and 37 GHz. 
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Fig. 9. — Radio light curves of PKS 0528+134. The vertical yellow shaded band corresponds 
to the interval when MDM optical and XMM-Newton X-ray observations of this source were 
performed. 
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frequencies from 7-rays through radio indicates no significant variabihty in the 7-ray regime, 
moderate variabihty in the X-rays {\AF/F\ ~ 50 %) and most radio frequencies {\AF/F\ < 
20 %) on time scales of ~ 1 - 2 weeks, and variabihty of up to AR < 1™^^ in the optical bands 
on time scales of several hours. Our data are not sampled densely enough for a meaningful 
analysis of time lags among different bands. 



4.2. Spectral Variability Analysis 

In order to test whether the optical flux variability discussed in the previous section is 
associated with spectral changes, we evaluated color indices B - R and V - R as a measure of 
spectral hardness. These were calculated for any pair of B (V) and R magnitudes measured 
quasi-simultaneously, i.e., during the same night, by MDM and SMARTS. If conditions were 
good enough to extract more than one high-quality (error in magnitude < 0.1) V or B and R 
band data point per night, magnitude measurements taken within < 1 hr of each other were 
used to calculate color indices. The resulting flux - spectral hardness correlations, indicated 
by the color- magnitude diagrams (color vs. R magnitude), are shown by the open symbols 
with dashed error bars in Figure [TOl The data clearly indicate color variability. Specifically, 
a weak softer-when-brighter trend seems to be present in both diagrams. 

A correlation analysis of the color-magnitude data sets yields a Pearson's correlation 
coefficient of r = —0.50, for the V - R color vs. R magnitude correlation, and r = —0.45 
for the B - R color. In both cases the obtained correlation coefficients indicate a weak 
negative correlation between color and magnitude, which confirm the redder-when-brighter 
trend. In order to assess the significance of these color correlations, we performed Monte- 
Carlo simulations of 10 million data sets with the same number of data points and the 
same spread in values as our data. The R magnitudes and color indices are assumed to be 
randomly distributed and not correlated to each other. For each simulated data set, the 
Pearson's correlation coefficient was evaluated in the same way as done for our observational 
data. This resulted in a chance probabilities of P{< r) = 5 x 10^^ and P{< r) = 2 x 10^'', 
respectively, of obtaining a correlation coefficient more negative than the ones resulting from 
our observational data. This seems to provide strong evidence for the presence of a redder- 
when-brighter trend in PKS 0528+134. 

Such a redder- when-brighter trend has been observed in other FSRQs, e.g., 3C 454.3 



( iRaiteri et al.ll2008l ). where it has been partially attributed to a contribution of emission 
lines from the BLR to the B band. In order to test whether this may also be the cause of 
the color variability described above, we utilize the line fluxes inferred from the spectrum 
shown in Figure [T] We note that the red-shifted CIV line falls within the B-band range. At 
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Fig. 10. — Color-magnitude diagrams for B - R (top panel) and V - R (bottom panel) vs. 
R-band magnitude. Open symbols with dashed error bars indicate the original data; filled 
symbols with solid error bars show the data data after correction of all magnitudes for the 
contributions from the CIII] and CIV emission lines. 
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A = 4740 A, the B-band filter has a transmission coefficient of approximately fx = 70% of 
its maximum value. Hence, using the bandwidth of the B filter of Au = 1.4 x 10^^ Hz, we 
find that the CIV line will make an effective contribution of F^^jy = Fqiv fx/ Av ^ 2.1 /i Jy 
to the B-band flux. An analogous calculation of the contribution from the CHI] line to the 
V and R bands yields -FJ^ct//] ^ 1.0 /iJy and F^^jjj^ ^ 1.6 /iJy. We corrected all B, V, 
and R magnitude values for these line contributions and re-evaluated the color magnitude 
correlations. The resulting points are shown as solid symbols in Figure [TOl This leads to an 
overall slight shift towards fainter (larger R) magnitudes and redder B - R colors. However, 
the overall color variability trend and its significance remain unaffected by this correction. 

We also analyzed the spectral variability in the X-ray regime. For the XMM-Newton 
data, we analyzed the variability of the hardness ratio F2-iokev/ Fo.2-2keV- Figure [TT] illus- 
trates that the hardness ratio is consistent with being constant over the four XMM-Newton 
observations performed as part of the core campaign. The RXTE X-ray energy index a was 
analyzed over a period of 150 days as plotted in the top panel of figure [121 A fit to a constant 
results in = 1.31 is found, indicating very moderate spectral variability. The comparison 
of the XMM-Newton spectra with the Suzaku spectrum from 2008 September/October (see 
§3.2.3p suggests that the X-ray spectral index remains stable on very long time scales, even 
throughout substantial (factor ~ 2 - 3) flux variations. 

We also analyzed the Fermi 7-ray spectral index variability over the extended campaign 
period. As shown in the bottom panel of figure [T21 no significant variability is found. The 
best fit found for the spectral index in this period of time is F = 2.9 ± 0.2, which is slightly 
softer than the mean spectral index ( F = 2.64 ± 0.06) ob served over the first eleven months 



of Fermi data (dashed horizontal line lAbdo et al.ll2010al ). 



4.3. Optical Polarization Variability Analysis 

The polarization variability analysis on PKS 0528+134 was performed using the R- 
band polarimetric observations at the 1.8 m Perkins telescope of Lowell observatory during 
a twelve day period from October 15 to 26, 2009, and two measurements from the University 
of Arizona / Steward Observatory Fermi support program from the same period. As shown in 
figure [T3l the Perkins polarimetry data indicate strong variability of the degree of polarization 
{xl = 20.2 for a fit to a constant) and the electric vector position angle (EVPA). The 
variability of the polarization parameters occurs on time scales of ~ 1 - 2 days and shows 
no obvious correlation with the fluxes in the optical or any other wavelength band. 

The significant degree of polarization in the optical (and radio) hints towards a sub- 
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Fig. 11. — Light curves in soft (0.2 - 2 keV) X-rays (top panel), hard (2 - 10 keV) X-rays 
(middle), and the hardness ratio. No significant spectral variability was found. 
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Fig. 12. — Top panel: X-ray (2 - 10 kcV) energy spectral index a vs. time. Bottom: Fermi 
7-ray photon spectral index (Fph = a + 1) vs. time. The dashed horizontal line corresponds 
to an index of 2.64, the mean index observed over the first eleven months of Fermi data. 
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Fig. 13. — Simultaneous photometric and polarimetric data for PKS 0528+134 on seven 
days during October 15 - 26, 2009. Top panel: R-band magnitude. Middle panel: Degree of 
polarization. Bottom panel: Electric vector position angle (EVPA). While there is significant 
variability in the R-band fiux, degree of polarization, and EVPA, we did not find a significant 
correlation between polarization and fiux states. 
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stantial synchrotron contribution to the emission at these wavelengths. We will discuss 
implications of our polarization results in Section [71 



Structure of the Parsec Scale Jet 



The quasar PKS 0528+134 is monitored monthly by the Boston University (BU) group 
with the Very Long Baseline Array (VLBA) at 43 GHz within a sample of bright 7-ray 
blazars [^. The source was also included in a 2-week campaign of observations of 12 7- 
ray blazars organized in 2009 October when 3 additional VLBA epochs at 43 GHz were 
obtained (VLBA project S2053). Figure [TH shows the total and polarized intensity images 
of the quasar in 2009 Autumn. The VLBA da t a wer e calibrated, imaged, and modeled in 
the same manner as discussed in IJorstad et al.l (120051 ). As we did for the optical data, the 
p olarization parameter s were evaluated taking into account the statistical bias as detailed 



m lWardle fc KronbergI (1l974[ l. Table |5] gives the parameters (flux, position, size, degree and 
position angle of polarization) of the main features seen in the radio jet during this period. 
Figure [15] shows the light curve of the VLBI core at 43 GHz of the quasar over the last 
three years as monitored by the BU group. According to Figure [15] the parsesc scale jet of 
PKS 0528+134 was in a quiescent state in 2009 Autumn. The core was moderately polarized, 
P ~ 2 %, with a stable position angle of polarization in a range between 42° and 50° that 
aligns within 10° with the jet direction as determined by the position angle of the brightest 
knot C2 with respect to the core. The closest knot to the core, CI, has the highest level 
of polarization, P ~ 10 %, with the position of polarization perpendicular to the local jet 
direction. Although we do not observe superluminal knots in 2009 Autumn, the quasar is 
known t o have a very high apparent speed in the jet, implying a high bulk Lorentz factor, 
r ~ 30 (Ijorstad et al1l2005h . 



6. Spectral Energy Distributions (SEDs) Modeling 

In figure [161 we present the SEDs of PKS 0528+134 from radio to 7-rays (blue filled 
circles) corresponding to the four XMM-Newton observations (September 8, 10, 11, 14). 
The optical data have been dereddened assuming Ay=2.782, As = 3.62 and Aji = 2.24, 
as given in the NASA/IPAC Extragalactic Databasqj according to the sky dust map given 



http://www.bu.edu/blazars 
http: / / nedwww.ipac.caltech.edu/ 
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Table 5. Parameters of Jet Components 



Epoch 


Knot 


S[Jy] 


R [mas] 


en 


a [mas] 


P[%] 


X[°] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


16 Sep 


AO 


0.842±0.085 


0.0 




0.048 


2.5±0.8 


47±7 




CI 


0.14±0.025 


0.165±0.025 


71.5±0.5 


0.107 


12.3±2.2 


-12±8 




C2 


0.35±0.050 


0.70±0.05 


64.3±0.5 


0.522 


<2.0 




14 Oct 


AO 


0.702±0.065 


0.0 




0.000 


1.4±0.6 


49±8 




CI 


0.14±0.025 


0.136±0.025 


82.4±0.5 


0.187 


8.0±2.5 


0±10 




C2 


0.29±0.050 


0.71±0.05 


61.8±1.0 


0.568 


5.4±1.8 


93±5 


16 Oct 


AO 


0.794±0.055 


0.0 




0.035 


1.0±0.5 


50±8 




CI 


0.25±0.025 


0.127±0.025 


68.7±0.5 


0.250 


8.2±2.0 


-3±10 




C2 


0.37±0.07 


0.72±0.05 


63.4±0.5 


0.510 


4.0±1.5 


91±8 


20 Oct 


AO 


0.826±0.065 


0.0 




0.023 


2.1±0.6 


42±5 




CI 


0.21±0.025 


0.143±0.025 


69.4±0.5 


0.240 


<5 






C2 


0.36±0.05 


0.72±0.05 


64.2±0.5 


0.490 


4.4±1.5 


91±6 


28 Nov 


AO 


0.794±0.045 


0.0 




0.035 


<1.0 






CI 


0.25±0.02 


0.137±0.025 


68.7±0.5 


0.250 


<4 






C2 


0.37±0.05 


0.72±0.05 


63.4±0.5 


0.510 


<2.4 





Note. — Columns: 1 - epoch of the observation; 2 - component designation; 3 - flux of component; 4 

- distance of component from the VLBI core; 5 - position angle of component with respect to the core; 
6 - diameter of component; 7 - degree of polarization of component; 8 - position angle of polarization of 
component 



1.5 1.0 0.5 0.0 -0.5 -1.0 
Milliarcseconds 



Fig. 14. — 43 GHz total (contours) and polarized [color scale) intensity im- 
ages of PKS 0528+134 during 2009 Autumn. The highest contour corresponds to 
5'peak=850 mJy/beam, while the yellow color indicates the highest polarized flux of ^'p^^^^j^ = 
30 mJy/beam, for a beam of 0.33 x 0.15 mas^ at PA=-10°. Total intensity contours corre- 
spond to 0.25, 0.5, 64 % of the peak. Sticks over the polarized intensity contours indicate 
the plane of polarization. The designation of components corresponds to Table [51 
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Fig. 15. — Light curve of the VLBI core at 43 GHz. The vertical lines indicate the period of 
intensified monitoring in 2009 Autumn. 
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by ISchlegel et al.l ( 119981 ) . We found in sections 14.11 and 14.21 that the Fermi 7-ray flux and 
spectral index of PKS 0528+134 did not show significant variability in the interval including 
the core campaign. Accordingly, the same 7-ray spectrum was used in the four SEDs. As 
in its flare states, in this quiescent state the SED of PKS 0528+134 is characterized by 
two peaks, a low-energy peak between the far infrared and optical spectral bands, and a 
high-energy peak at MeV - GeV energies. As can be seen, in all the SEDs the high-energy 
component dominates the bolometric output by a large amount. 

We produce model fi ts to all four SEDs u s ing th e equilibrium version of the leptonic one- 



zone model de yeloped bvlBottcher fc Chiang (120021 ). This equilibrium model is described in 



more detail in lAcciari et al.l ( 120091 ). and we here summarize its main features. The observed 
electromagnetic radiation is interpreted as originating from ultrarelativistic electrons (and 
positrons) in a spherical emission region of co-moving radius R, which is moving with a 
relativistic speed /3rc, corresponding to the bulk Lorentz factor F. Depending on the viewing 
angle 6 between the jet direction and the line of sight, the transformations of photon energies 
and fluxes is characterized by the Doppler factor D = (F[l — /3rcos^])~^. The size of the 
emission region is constrained by the shortest observed variability time scale (5tvar,min through 

R < C5tvar,min^/(l + z) < 8 X lO^^ ((5tvar,min/d) (D/IO) CUi. 

Ultrarelativistic electrons are assumed to be instantaneously accelerated at a height zq 
above the accretion disk into a power-law distribution in electron energy. Eg = jmeC^, at a 
rate per unit volume and unit Lorentz factor interval given by ^(7) = <5o7~'^ with a low- and 
high-energy cutoffs 71 and 72, respectively, and injection spectral index q. An equilibrium 
between this particle injection, radiative cooling, and escape of particles from the emission 
region yields a temporary equilibrium state described by a broken power-law. The time 
scale for particle escape is parameterized through an escape time scale parameter r^esc > 1 
as tesc = VescR/c- The balance between escape and radiative cooling will lead to a break 
in the equilibrium particle distribution at a break Lorentz factor %, where tesc = ^cooi(7)- 
The cooling time scale tcooi is evaluated self-consistently taking into account synchrotron, 
synchrotron-self- Compton (SSC) and external Compton (EC) cooling. The number density 
of injected particles is normalized to the resulting power in ultrarelativistic electrons 
propagating along the jet. The magnetic field B in the emission region is pre-specified as a 
free parameter. It corresponds to a Poynting flux along the jet. Lb = nR^ T^(3r cu'^ where 
u'b = -B^/(87r) is the magnetic field energy density in the co-moving frame. For each model 
calculation, the resulting equipartition parameter, 63 = Ls/Lg is evaluated. 

Once the quasi-equilibrium particle distribution in the emission region is calculated, our 
code evaluates the radiative output from synchrotron emission, SSC, and EC emission self- 
consistently with the radiative cooling rates. If the occasional indication of a blue bump in 
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the optical spectrum can be associated with the accretion disk, we can estimate an accretion 
disk luminosity from the corresponding approximate uF^, flux in the UV regime of uF^^^^^ ~ 
8 X 10^*^ Jy Hz as L/) ~ 1.7 x 10^^ erg s~^, which we use for our model flt s. The accretion disk 
emissi on is modelled as a multi-color blackbody spectrum according to a lShakura fc Sunyaev 
fll973h disk model. 

In addition to direct accretion disk emission, external radiation may originate as line 
emission from the Broad Line Region (BLR). We can est imate the tot a l lurn inosity of the 
BLR line emission using the bright quasar template of iFrancis et al.l (119911 ). normalized 
to the observed value of the CIV emission line luminosity of Lciv = 2.5 x 10^^ erg s~^. 
Substantial contributions to the BLR luminosity, in addition to GUI] and CIV observed here, 
are expected to arise from Fe II, Lya and Ly/3, H/3 and H7, Mg II, and He II, among others. 
The total BLR luminosity is expected to be Lblr ~ 7.1 Lciv ~ 1-8 x 10^^ erg s^^. This is an 
order of magnitude lower than the value we adopt for the accretion disk luminosity. However, 
the relevance to external Compton scattering depends on the photon fleld energy density in 
the rest-frame of the emission region and hence on the geometry and the Doppler boosting 
of external photons into the co-moving frame of the emission region. The energy density of 
the direct disk emission in the AGN frame is u^isk = Ljj/^inzlc) where Zq is the distance 
of the emission region from the central supermassive black hole. If zo 3> Rd T^, where 
Rd is the characteristic radius of the annulus of maximum energy output of the accretion 
disk, accretion disk photons will be strongly red-shifted in the emission-region rest frame. 
However, in the near-fleld regime, zq -C Rd T^, even the accretion disk photons entering 
the emission region from behind in the AGN rest frame will still be blue-shifted and their 
energy density enhanced in the emission region rest frame (see, e.g., iDermer fc Schlickeiser 



I993I ). If the emission region is located within the characteristic radius of the BLR, Rblr, 



the BLR radiation fleld can be treated as approximately isotropic, and it will be blue-shifted 
into the emission-region rest frame, and enhan ced by a factor of ^ compared to its AGN 
rest-frame value of ublr = Lblr/ {4:71 R^^j^j^c) ( jSikora et al.lll994l ). 

Lacking knowledge of the size of the BLR and the precise location of the emission 
region in PKS 0528+134, both direct disk and BLR emission contributions remain plausible 
as dominant sources of external photons for the external- Gompton process. In the following, 
we have chosen our model parameters (in particular, zq) such that we expect the direct disk 
emission to dominate. As we will see below, this allows for satisfactory flts to the SED and, 
in particular, the Fermi-LAT 7-ray emission from PKS 0528+134. Therefore, in order to 
avoid the introduction of another unconstrained parameter, Rblr, we restrict our modeling 
efforts to using the direct accretion disk emission as the dominant contributor to the external 
radiation fleld to evaluate the external radiation Gompton (ERG) component of the SED. 
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Frequency Hz 



Fig. 16. — Spectral energy distribution (SEDs) of PKS 0528+134 corresponding the first 
XMM-Newton observation. Blue points correspond to the observational data. The radio 
data correspond to the integrated radio flux from single-dish measurements. The orange 
radio point at 43 GHz is the flux of the VLB A core component. The continuous lines 
correspond to the best flt for each SED. The dotted line represents the synchrotron spectrum, 
the dot-dot-dashed line corresponds to the Synchrotron self Compton (SSC), the dot-dashed 
line is the external radiation Compton (ERG), and the short-dashed line is the disk thermal 
component. 
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As a representative example, the fit to the SED of PKS 0528+134 during the first 
XMM-Newton observation is shown in Figure dH The parameters used for this particular 
fit are shown in table O The synchrotron (dotted line), synchrotron self Compton (SSC) 
(dot-dot-dashed line), external radiation Compton (ERC) (dot-dashed line), and the direct 
disk emission (short-dashed line) components are shown separately, in addition to the total 
SED fit curve. 

We were able to achieve good fits for all four SEDs of PKS 0528+134 in quiescence. 
Table [7] lists the most relevant parameters used in each fit. Fit parameters close to equipar- 
tition could be found for all four SEDs. In general, no obvious correlation between the 
different parameters was found. However, a strong correlation (Pearson's correlation coeffi- 
cient r ^ 1) was found between the magnetic field and the optical flux (R-band), reflecting 
the synchrotron dominance in the optical band. 

We point out that our model only includes the emission from the blazar zone, assumed 
here to be on sub-pc scales. It is expected that the radio emission originates in the more 
extended (pc to kpc scales) jet which is not included in our model. Therefore, our fits under- 
produce the radio spectra in all SEDs. While most of the radio data shown in Figure [11] 
are obtained by single-dish instruments and therefore represent the integrated flux over all 
radio components, we have also included the 43 GHz flux from the VLBA core component 
on 16 September 2009, just 2 days after the last XMM-Newton observation. Even this core 
radio flux is under-represented by our model, suggesting that the higher-frequency emission 
originates on even smaller scales than the radio core. 



Discussion 



As mentioned in the previous section, good flts with a one-zone leptonic SSC + ERC jet 
model were possible with parameters close to equipartition. Our flt parameters, in particular, 
the magnetic flelds of i? ~ 2 - 3 G, characteristic electron energies, 7mi n ^ 10^, and jet 
power s of L, 



]^q45 gj,g g 1^ g^j^g rough agreement with the flt results of iMukherjee et al. 



( 1999 ). However, we need to point out that the model used in Mukherjee et al. ( 19991 ) is not 
precisely the same as used in this paper , as it was based on a time-average of an evolving 
particle distribution along the jet. Also, iMukherjee et al.l ( 1l999l ) used a different cosmology 
(a matter-dominated Universe with go = 0.5), resulting in a substantially different luminosity 
distance. Therefore, our results are not directly comparable. We notice that the Fermi LAT 
7-ray spectra of PKS 0528+134 are systematically softer th an the photon indices Pph ~ 2.2 
- 2.6 found during the EGRET era (IMukherjee et al.lll999l ). Therefore, our flts to the LAT 
SEDs require signiflcantly steeper particle spectral indices. 
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Table 6. Parameters used in the fit for tlie SED of PKS 0528+134 corresponding to the 
XMM-Newton observation of September 8, 2009 (JD 2455082.94). 



Parameter 


Value 




1400 


'ymax 


10^ 


Injection electron spectral index 


3.65 


Escape time parameter (fcsc = '/osc R/c) 


)7csc = 50 


Magnetic field [G] 


2.05 


Injection height [pc] 


ZQ = 0.13 


Bulk Lorentz factor V 


20.4 


Accretion disk luminosity [10^^ erg/s] 


Ue = 17 


Blob radius [cm] 


1.9 X 10^6 


Black hole mass [Mq] 


1.08 X 10^ 


Observing angle [degrees] 


e^bs = 3.02 


Doppler factor 


D = 19.0 


Redshift 


Z = 2.06 


Le (jet) [erg/s] 


2.3 X 10"'^ 


Lb (jet) [erg/s] 


2.3 X 10''^ 




1.0 



Table 7. Relevant fit parameters for the SEDs of PKS 0528+134 



SED 


B [G] 


Fr [Jy Hz] 


Fx [erg cm s ^] 


r 


q 


»? 




Lb 


Lb /Le 


20 [pc] 


1 


2.05 


1.8399e+ll 


1.6249e-12 


20.4 


3.65 


50 


2.30e+45 


2.30e+45 


1.0 


0.13 


2 


2.06 


1.6570e+ll 


1.3224e-12 


20.1 


3.85 


50 


2.39e+45 


2.26e+45 


0.95 


0.13 


3 


3.7 


3.8000e+ll 


1.3702e-12 


17.5 


3.7 


70 


1.18e+45 


1.73e+45 


1.46 


0.11 


4 


2.21 


1.9799e+ll 


1.4367e-12 


21.1 


3.65 


40.3 


2.24e+45 


2.40e+45 


1.07 


0.135 
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We also need to caution that the model contains a large number of poorly constrained 
parameters, and in many cases, different parameter combinations might be able to produce 
similarly acceptable fits. Therefore, conclusions about correlations between model parame- 
ters and observables may not be unique. 

A simi l ar stud y of 7-ray bright blazars in quiescence has recently been published by 



Abdo et aLlteOlObl ). Those authors observed 5 blazars (PKS 0208-512, Q 0827+243, PKS 1127- 
145, PKS1510-089, and 3C 454.3) in their low-activity state with Suzaku and Swift in X-rays 
and optical/UV, and analyzed the simultaneous Fermi-LAT data. All of those blazars showed 
X-ray continua consistent with a hard (Fph ~ 1.5) single power-law spectrum, in agreement 
with our XMM-Newton and Suzaku results on PKS 0528+134. The broadband SEDs of all 
five blazars were dominated by their high-energy (7-ray) output, as in PKS 0528+134. In 
contrast to our XM M-Newton results on PKS 0528+134, three of the five blazars observed by 



Abdo et all (l2010bl ) (PKS 0208-512, PKS 1127-145, and PKS 1510-089) did show significant 



short-term vari ability on time scale s of ~ 5 - 10 hr. However, we need to point out that the 



observations by I Abdo et al.l ( 1201 Obi ) did not adhere to the strict quiescence criterion required 
for triggering our PKS 0528+134 campaign (i.e., the Fermi-LAT 7-ray flux remaining below 
the lowest EGRET flux or upper limit persistently for at least two weeks, see §T]). In fact, 
some of their targets (in particular, PKS 0208-512, PKS 1510-089, and 3C 454.3) exhibited 



substantial 7-ray flux and even 7-ray flares within just a 



Therefore, the X-ray variability reported in I Abdo et al.l (l2010bl ) may not be characteristic of 



'ew day s of the Suzaku observations. 



a truly quiescent state of those blazars as investigated in the case of PKS 0528+134 reported 
in this paper. 

In our analysis of optical spectral variability we found a weak anti-correlation between 
the B-R color and the R-band magnitude. A possible explanation for this softer when brighter 
trend can be found in the interplay between the two radiation components that contribute 
photons to the optical flux: the synchrotron component, which is generated in the jet itself, 
with a steep spectrum, always dominating at low (R-band) frequencies, and slowly varying 
emission components associated with the accretion disk and the BLR. The direct accretion 
disk emission is expected to peak in the ultraviolet and may therefore contribute at the blue 
end of the optical spectrum. In Section |321 we have shown that the observed color variability 
can not be attributed to the contribution from the CIII] and CIV emission lines in the optical 
spectrum of PKS 0528+134. Therefore, the obvious candidate is direct, thermal accretion 
disk emission. The suggestive trend of increasing polarization with increasing wavelength 
lends further support to this hypothesis. 

Given the evidence we found for a substantial accretion disk with a luminosity of ~ 
1.7 X 10'^^ erg s~^, it is worth investigating whether at least part of the observed X-ray 
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emission may result from Comptonization of soft disk photons in a hot, thermal corona above 
the accretion disk. In the extreme scenario, in which the entire XMM-Newton spectrum is 
produced by the corona. Figure Hn] illustrates that one would require the energy dissipated 
in the corona to be of the same order as that dissipated in the optically thick disk. In 
the case of Comptonization of soft photons with energy far below the X-ray regime, the 
X-ray spectral index Fph ~ 1.6 can then be translated into a Comptonization parameter 
y = 16 Q^tt {I + tt), where 9 = kT/^m^c^) is the dimensionless coronal temperature, 
through 



yielding a value oi y = 1.85. This would imply coronal parameters of tt (1 + tt) = 0.116. 
However, in order for the X-ray spectrum to be dominated by coronal Compton emission, 
the jet SSC emission would have to be strongly suppressed compared to the model fits 
we presented here. This could, in principle, be achieved by the choice of a much larger 
emission region size R. Such a choice might be problematic for two reasons: (1) For R > 
1.7 X 10^^ {D/20) cm, causality would not allow for variability on time scales of a day or less; 
(2) as the observed synchrotron emission would require the same magnetic field as chosen in 
our SED fits, this would result in the energy density in the leptonic particle population being 
far below equipartition with the magnetic field. While both arguments may not strictly rule 
out such a scenario, we strongly prefer our fit scenario, close to equipartition, and allowing 
for ~ day scale variability. 

We found that the shortest observed variability occurs at optical wavelengths. In the 
standard interpretation, this is the region near the high-frequency end of the synchrotron 
spectrum, emitted by the highest-energy electrons with the shortest radiative cooling time 
scales. The synchrotron cooling time scale for electrons radiating via synchrotron emission 
in the V band, is 



Using the redshift of z = 2.07, a characteristic magnetic field of i? ~ 2 G, and a typical 
Doppler factor of D = 19 from our fit results, we find an observed synchrotron cooling time 
of r°y^'(V") ^ 0.9 hr. Since the SED of PKS 0528+134 is strongly dominated by the 7-ray 
emission, the actual radiative cooling rate is expected to be shorter than the rate expected 
from synchrotron emission alone, by a factor corresponding to the Compton dominance (the 
ratio of power output in the high-energy vs. the synchrotron component), which is of order 
10. Therefore, the total radiative cooling time scale of electrons emitting synchrotron in the 




(2) 



T°^\V)^2.2 X 10^(l + z) 



(3) 
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V band is likely of the order of ~ 1 minute or shorter. This suggests that the light curves 
of variability in all optical bands are strongly dominated by light-travel time effects rather 
than microphysical (acceleration and cooling) time scales. 

The very moderate degree of variability observed in all wavelength bands, in tandem 
with the short radiative cooling time scales, implies on-going particle acceleration in the 
quiescent state of PKS 0528+134. On the other hand, the significant polarization variability 
hints towards a turbulent process. Steady, on-going particle acceleration can be envisioned 
in shear-layers of radially stratified jets or particle acceleration at standing features, such as 
re-coUimation shocks. However, the magnetic fields in shear-layers of radially stratified jets 
are expected to be highly ordered due to the directionality of the differential motion in those 
systems. Therefore, our results seem to favor scenarios involving turbulent acceleration, 
possibly associated with stationary re-collimation shocks. 

Further evidence for turbulent processes comes from our polarization results. A power- 
law distribution of non-thermal electrons with a spectral index p can produce synchrotron 
emission with a maxin ium polarization degree of P^'^^ = {p + l)/{p + 7/3) for a perfectly 



ordered magnetic field (IRybicki fc Lightmanlll986l ). Based on our modeling results in Section 



[6l the high-energy (cooled) part of the electron spectrum is expected to have a spectral index 
of p = g + 1 ~ 4.6. This would yield a maximum polarization degree of P™^^ ^ 70 %. The 
fact that the actual degree of polarization remains below ~ 10 % indicates the magnetic 
field is tangled on size scales much smaller than the size of the emission region. The rapid 
(day scale), apparently random variation of the degree of polarization and the EVPA seems 
to indicate a turbulent process resulting in a large n umber of individual cells with randomly 



oriented magnetic fields. Monte Carlo simulations (jD'Arcangeld l2010l ) indicate that for a 
100 % chaotic magnetic field (i.e., no ordered field component), 150 turbulence cells result 
in a degree of polarization of p = (6 ± 3) %, in agreement with the observed values for 
PKS 0528+134. Given a typical variability time scale in the optical of ~ 1 day, the size 
of individual turbulence cells can be estimated as i?ceii ~ cSt^g.^/ {N^^^ D) ~ 2 x 10^^ cm 
for a Doppler factor of Z) ~ 20 (see ^ and = 150 turbulence cells. Comparing the 
optical polarization variability to the polarization variability of the radio core (see one 
can speculate that when the optical polarization is high (> 5 %), the optical EVPA seems 
to be more closely aligned with the EVPA of the radio component closest to the radio core. 



8. Summary and Conclusions 



Over the last two decades PKS 0528+134 has become an important target for multi- 
wavelength observations because of its high luminosity from radio through 7-rays and the 
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extreme flux and spectral variability that it shows in its flaring states. In this paper, we 
have presented multiwavelength observations of PKS 0528+134 involving the XMM-Newton, 
RXTE, Suzaku, and Fermi satellites as well as many ground based radio and optical tele- 
scopes. Our main goal was to characterize this 7-ray loud quasar in a quiescent state to 
improve the understanding of SEDs and variability patterns of this source, and of blazars in 
general. The variability analysis of the collected data, and the construction and modeling of 
four SEDs of this source, yielded the following results: 

• No signiflcant short term flux and spectral variability (as determined by data in the core 
campaign) was found in 7-rays, X-rays and most radio bands. However, for the same 
time interval, signiflcant flux variability with AR < 1™^^ on time scales of several hours 
was found in the optical, accompanied by a weak spectral softening with increasing 
flux. The latter trend may be interpreted as a steady contribution of the accretion disk 
flux at the blue end of the optical spectrum. Optical spectropolarimctry suggests an 
increasing degree of polarization towards longer wavelengths, lending further support 
to the hypothesis of synchrotron emission contributing an increasing fraction of the 
emission towards the red end of the spectrum. 

• Data analysis based on a more extended interval (two months or more) shows no 
significant 7-ray fiux and spectral index variations, but moderate flux variability in 
the X-rays {\AF/F\ ~ 50 %) and at radio {\AF/F\ < 20 %) frequencies on time scales 
of ~ 1 - 2 weeks. A comparison between XMM-Newton spectra of 2009 September and 
the Suzaku spectium of 2008 September/October suggests that the X-ray spectral index 
remains stable in the quiescent state of PKS 0528-1-134 even throughout substantial 
(factor ~ 2 - 3) long-term flux variations. 

• We constructed four SEDs of PKS 0528+134 in quiescence. Our results show that 
even in the quiescent state, the bolometric luminosity of PKS 0528+134 is strongly 
dominated by its 7-ray emission, although the 7-ray spectra are signiflcantly steeper 
than found during the EGRET era. 

• We fitted the four SEDs with a leptonic combined SSC+ERC jet model. In this 
model, the low energy component is produced by the sum of the synchrotron process 
in the jet and the disk luminosity, and the high energy emission is due to the sum of 
the synchrotron self-Compton and the external radiation Compton contributions. Fit 
parameters close to equipartition were found for all SEDs. 

• The moderate variability in most wavelength bands, compared to the expected short 
radiative cooling time scale, implies the persistence of particle acceleration on long 
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time scales. This may favor acceleration scenarios based on standing features, such as 
re-coUimation shocks. 
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